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ABSTRACT 

The line ratios from coronal gas in Seyferts can be successfully fitted with pho- 
toionized clouds of high densities and low volume filling factor. The ionization pa- 
rameter implied is sufficiently high that models must consider the effect of radiation 
pressure from the active nucleus. In spite of the nucleus gravitational force, radia- 
tion pressure is sufficiently strong to compress and radially accelerate the internally 
stratified gas clouds provided these contain small amounts of dust (~ 10% of so- 
lar neighborhood value). This radial acceleration could explain the blueshift of the 
coronal lines relative to systemic velocity without the need of invoking an ambient 
'pushing' wind. Embedded dust has the interesting effect of making the photoionized 
clouds marginally ionization bounded instead of matter bounded. 
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1 INTRODUCTION 

Studies of the emission lines in Seyfert galaxies have indi- 
cated the presence of two different regions: a region of very 
high density with large velocity spread - the broad line re- 
gion (BLR), and a region of low gas density with a much 
smaller velocity spread - the narrow line region (NLR). 
However, there is also a set of forbidden lines of extremely 
high excitation, the so-called 'coronal' lines which in- 
cludes for instance [Fevn] A6086, [Fex] A6374, [S vm] A9913, 
[SiVl] 1.963/xm and [Sivn] 2.483/mi. Most models favor pho- 
toionization over collisional ionization as excitation mecha- 
nism (e.g. Korista & Ferland 1989; Oliva et al. 1994; Moor- 
wood et al. 1996). Earlier models considered that the filling 
factor of the coronal gas approached unity with a density 
lower than that of the NLR. The advent of far infrared mea- 
surements of the density sensitive [Nev] 14.3 ^im/24.3 )im 
line ratio with ISO (Moorwood et al. 1996: Mo96), however, 
suggests much higher densities of ~ 5 000 cm -3 implying a 
very low volume filling factor for the coronal gas. In Paper I, 
Binette et al. (1997) proposed that these lines originate from 
individual gas clouds - as for the NLR - but characterized by 
an unusually high ionization parameter, Uo [cf. equation (1) 
below]. A similar gas geometry underlies the extensive grid 
of photoionization models of Ferguson, Korista & Ferland 
(1997). In Paper I, the clouds responsible for the coronal 
lines are always matter-bounded while they are essentially 
ionization-bounded in the models of Ferguson et al. (1997) 
or in the photoevaporing dusty cloud model of Pier & Voit 
(1995). One advantage of the matter-bounded clouds is that 
they are easier to accelerate by radiation pressure owing 
to their lower mass. Since the coronal lines are known to be 



systematically blueshifted in Seyfert galaxies (Penston et al. 
1984) relative to the systemic velocity (by 35kms _1 in the 
case of Circinus, see Oliva et al. 1994; hereafter OSMM), 
our justification of using matter-bouneded clouds was to ac- 
count for this property by incorporating self-consistently the 
effect of radiation pressure in our cloud model. 

The ionization parameters required to fit the very high 
excitation coronal lines is very high, UoiJ 0.2, implying that 
the pressure exerted by the ionizing radiation exceeds the 
gas pressure (see Paper I). This has two interesting conse- 
quences: first, the photoionized cloud can be radially acceler- 
ated until the ram presure it exerts on the ambient medium 
equals the pressure of the latter, this would account for the 
blueshift of the coronal lines for certain geometries, second, 
a density gradient arises within the clouds due to radiation 
pressure (a distributed force) and as shown in Paper I, the 
density stratification generated in this way leads to a wider 
range in excitation of the ionization species present in the 
cloud. The question addressed in this Letter is whether the 
pressure from the nuclear radiation field is sufficient to cause 
a strong radial acceleration of the clouds despite the gravi- 
tational pull exerted by the stars and the nuclear black hole. 



2 THE MODEL 

The multipurpose code mappings IC was used to compute 
the photoionization models. The atomic data is taken from 
a compilation of Ralph Sutherland (cf. Appendix in Ferruit 
et al. 1997) . The uncertainties in collision strengths of many 
coronal lines are often large (see review by Oliva 1997 and 
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the discussion of the 'iron conundrum' by FKF) and may 
dominate the errors in the predicted line strengths. 

2.1 Gas abundances and dust content /id 

As is customary and in the absence of other information, 
the abundances adopted will be solar (Anders & Grevesse 
1989) but depleted according to the dust content. The lat- 
ter is defined by the quantity /id which is the dust-to-gas 
ratio of the plasma expressed in units of the solar neigh- 
borhood dust-to-gas ratio. The effects of dust on the ther- 
mal and ionization structure are taken into account (Binette 
et al. 1993). However, a different depletion scheme is defined 
whereby the destruction of dust grains is assumed to return 
uniformly the depleted elements ('uniform return'). For in- 
stance, supposing /iD=0.1 (i.e. 10% of the solar neighbor- 
hood dust content), the gas phase abundances are derived 
from the fully depleted set of abundances (/id = 1) to which 
we restitute 90% of what has been depleted. This means 
that even for the heavily depleted elements like Ca, the gas 
phase abundance when /xd=0.1 is 90% solar. The justifica- 
tion behind this is that small values of fio are not the result 
of a different dust formation history but rather that of re- 
cent destruction (of normally formed dust) by sputtering, 
sublimation or optical erosion. 

2.2 The ionizing continuum 

We adopt a similar ionizing continuum as in Paper I, namely 
a powerlaw of index 7 = —1.3 (</?„ oc u +1 ) for the infrared- 
UV domain. This powerlaw joins at 2000 eV with a flatter 
powerlaw of index 7 = —0.7 to cover the X-ray domain. 
We impose a high energy cut-off of 100 keV(see Mathews 
& Ferland 1987). As is customary, we define the ionization 
parameter Uo as the ratio between the density of impinging 
ionizing photons and the gas density at the face of the cloud 

Uo = — = -52- , (1) 

cno cno 

where ip v is the monochromatic energy flux impinging on the 
cloud, v\ the Lyman frequency, qo the number of ionizing 
photons incident on the slab per cm 2 per second, c the speed 
of light and no the total gas density at the irradiated surface 
of the slab (all quantities relevant to the irradiated surface 
carry as a superscript or subscript). 

2.3 Acceleration and pressure stratification 

We adopt the formalism developed by Mathews (1986) and 
consider that our radially accelerated cloud has achieved in- 
ternal hydrostatic equilibrium in the radial (outward) direc- 
tion in the noninertial accelerating frame. Therefore, across 
a cloud (hereafter approximated as a slab) of total column 
density Nt and situated at a distance r from the black hole 
(BH), we have 

1 dP 

Iat = 9rad(N) - {g puU + a c i) = g ra d{N) - a , (2) 

where mn is the proton mass, g ra d(N) the local radiative ac- 
celeration due to radiation pressure, g pu li the gravitational 
acceleration due to the central massive BH and the inner 



stars, and a c i the instantaneous acceleration of the entire 
cloud. The volume force exerted by radiation pressure is 
simply F ra d — nmHg r ad(N) with its numerical expression 
in terms of photoelectric cross sections defined in Appendix 
A. 3 of Paper I. The clouds geometrical thickness is consid- 
ered negligible as compared to nuclear distance r, a good 
approximation in the current context of moderately thick 
clouds. This implies that for a cloud at r, the modified ac- 
celeration a' (= g vu u + a c i) is constant across its thickness 
Nt. We can regard a' as an eigenvalue for the acceleration 
of the cloud. In our scheme, the systematic shiftQ towards 
the blue of the lines requires that a' > g pu ii, that is a c j > 
in equation (2). The blueshift can either be the result of an 
earlier acceleration of the clouds (the clouds have reached 
their terminal velocity against the opposing drag and grav- 
itational forces), or of ongoing cloud acceleration. 

2.4 Implementation in mappings IC 

Independently of whether the acceleration corresponds to 
bulk motion or to inertial acceleration, a' is an important 
parameter for the cloud's density structure. In our model, we 
express this input parameter in terms of the dimensionless 
quantity $0 defined as ^0 = a! '/ 1 g® a d with g® ad the radiative 
acceleration at the irradiated face. 

As described in Appendix of Paper I, MAPPINGS in- 
tegrate equation (2) as well as simultaneously solve for the 
local ionization and thermal balance as a function of depth 
N. If Uo is sufficiently high, equation (2) implies a posi- 
tive pressure gradient within the slab therefore leading to 
a strong density gradient across the photoionized slab. We 
define the relative pressure difference between the irradiated 
face and the back of the slab, SU, as 

5H = (P b a :f - P 9 ° as ) /MIN (P b 3 :f , P 9 ° as ) , (3) 

where the denominator represents the isotropic pressure, P, 
of the ambient cloud-confining medium. Note that SU is 
approximately the ratio of dynamical pressure to ambient 
pressure. If SU < 0, the cloud is experiencing a dynamical 
pressure at the inner irradiated surface ('pushed' clouds by 
an ambient wind), while if ST1 > the clouds are suffering 
a drag force at their outer back surface as they are 'push- 
ing' against the ambient medium. All the models presented 
below correspond to this latter case of SU > 0. 

For the dust-free matter bounded slab discussed in Pa- 
per I where we neglected inertial and bulk acceleration 
(^0 = 0), radiative pressure was shown to induce a pressure 
difference in the slab as high as SI1 — 1.6 when Uo=0.5 . 
Since any realistic nuclear environment will imply a nonzero 
a' as a result of the gravitational field, it is the object of this 
study to investigate the effect of a positive ^o- 

In the general case of "l/o > 0, because g ra d decreases 
monotonically with depth as the impinging radiation is grad- 
ually absorbed, the right term of equation (2) become neg- 
ative at very large depths. This can result in <5n < de- 
pending on the choice of Nt, /xd and ^o- Although we could 
interpret this as being caused by ram pressure of an out- 
flowing ambient medium ('pushed' cloud), this domain of 

* It is believed that we see one hemisphere only (i.e. one ioniza- 
tion cone) of the NLR in both Circinus and NGC 1068. 
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parameters has not been explored: first, because the aim of 
the Letter is to relate the line blueshift to the acceleration 
caused by radiation pressure and not to an extraneous alter- 
native mechanism, second, because all the models considered 
below are not sufficiently thick to obtain <5II<0. 

To maximize the acceleration a', it is useful to increase 
$0 as much as possible. However, it was found that $o=0.5 
which is adopted hereafter cannot be exceeded significantly 
(otherwise 511 shrinks considerably). In effect, 1 3/o=0.5 main- 
tains two essential effects: it still allows a significant density 
gradient (therefore optimizing the line ratio fit, cf. Paper I) 
and it results at the same time in a strong acceleration a' . 
Clouds submitted to significant radiation pressure are par- 
ticularly prone to derimming through the establishment of 
a lateral flow as shown by Mathews (1986). For a review of 
the problem of cloud stability under conditions appropriate 
to the NLR, see Mathews & Veilleux (1989). 



3 THE EXCITATION OF THE CORONAL GAS 
IN CIRCINUS AND NGC 1068 

All detailed studies of the coronal gas reveal that the emit- 
ting gas encompasses many stages of excitation. This indi- 
cates that this gas must be significantly thick to the ionizing 
radiation rather than very optically thin. We can envisage 
two geometries, one broadly spherical and filled with low 
density gas: the onion ring geometry (Korista & Ferland 
1989; OSMM) the other consists of an ensemble of clouds 
of very low filling factor and therefore similar to the NLR 
gas in general except for its much higher excitation. What 
is clear is that although the range of ionization species of 
the coronal gas is wide, it does not encompass the very low 
excitation species of [On], [Oi], [Sll]... since the latter ap- 
pear to be kinematically distinct as shown by OSMM and 
Marconi et al. (1996). In Paper I, this was interpreted as an 
indication that the coronal line emitting clouds were mat- 
ter bounded as was also thought to be the cas e for the high 



excit ation clouds of the ENLR and NLR (?; Wilson et al 



1997) 



The thickness of the clouds is not defined arbitrarily but 
can be empirically set by the ratio between the highest and 
the lowest excitation species belonging to the coronal gas. 
In Paper I it was found that the lines of [Sivi] 1.963/^m and 
[Siix] 3.935pm worked very well for this purpose. The same 
procedure is adopted here whereby the slab is truncated at 
the point where the modeled [Si vi]/[Si ix] ratio matches the 
observed value. It is interesting to note that in the dust-free 
case, despite the difference in excitation between the two 
Seyferts, this truncation always seems to occur after ~ 60- 
65% of the ionizing photons have been absorbed. 

(i) The Circinus galaxy. 

If we repeat the dust-free calculations of Paper I with 
Uo=0.5 and density no=1200 cm -3 but set $0=0.5, we find 
that the modified acceleration is only a'= 2.9 x 10~ 7 cms -2 



- is twice that value) . As shown later in Section 



4.1 



(i-e- grad 

this is at least an order of magnitude below the acceleration 
due to the gravity exerted by the stars within 100 pc of the 
nucleus of the Milky Way. It would therefore appear that 
radiation pressure is insufficient in accelerating radially the 
line emitting gas. However, if we allow trace quantities of 
dust to be present in the cloud, radiation pressure is much 



a 



Mgv. 



U = 0.50 



i ■ 

0.9 - 
0.8 - 
15' 0.7 - 

£0.6 " 
\ 0.5 - 



-a 0.4 - 

| 0.3- 
O 




Model A : Circinus Galaxy 

— I 1 1 I 



_ VIII 

Cao 



U = 0.125 



Model B: 
NGC 1068 



i - 

0.9 - 
0.8 - 
0.7 ■ 
0.6 - 
0.5 - 
0.4 - 



Fe 



Nei 



50 100 150 200 350 300 

Ionization Energy (eV) 



Figure 1. Observed divided by modeled line flux ratios for (a) the 
Circinus galaxy and (b) NGC 1068. All ratios are expressed rela- 
tive to [Si ix] 3.935/Ltm. Both models A and B have no=1200 cm -3 
and a dust content of 10% of the solar neighborhood dust-to- 
gas ratio (fi£>=0.l, see text and Table. 1). The lines plotted in 
(a) include [O iv] 25.90Atm, [Si Vl] 1.963/wn, [Sivn] 2.483/im, 
[Siix]3.935/Ltm, [S rv] 10.54/Ltm, [S vm] A9913, 
[S IX] 1.25jtm, [NeV]14.32Atm, [Ne v] 24.31/mi, 
[Fe vn] A6086, [Fe x] A6374, [Fe xi] A7892, 
[Mgvn] 5.51/im and [Mgvm] 3.03/xm. 



[Cavm] 2.32yum, 
[Nevi] 7.66/xm, 
[Mgv] 5.62/mi, 



higher and dust absorption even dominates the photoelectric 
absorption terms contained in g ra d- For instance, if we set 
/iD=0.1, we obtain a'— 2.6 x 10 -6 cms -2 . The modeling of 
the observed line ratios reported by Mo96 and OSMM is 
presented in Fig. la (Model A) using fir>~0.1 . The fit is of 
comparable quality to that presented in Paper I or in Mo96. 
As can be expected, the effect of 'uniform return' depl etio n 
on the [Cavm] point is less than 0.05 dex (see Section 2.1). 
The resultant pressure stratification is somewhat smaller at 
<5n=1.2 . It was found that arbitrarily increasing either /xd 
or $0 to higher values degraded significantly the fit to the 
Circinus line ratios. 

The hypothesis of internal dust is also present in the com- 
prehensive dusty cloud model of Pier & Voit (1995). In their 
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model, the strong radiation pressure exerted on the dusty 
gas (situated near the inner molecular torus) leads to pro- 
gressive cloud evaporation and to the creation of an X-ray- 
heated wind. The smaller dust content and cloud's thick- 
ness of the matter-bounded clouds presented in this Letter 
leads to a smaller extinction (Ay— 0.15) than in Pier & Voit 
(1995). 

(ii) The Seyfert galaxy NGC 1068. 

We follow a similar procedure in the case of NGC 1068. 
The observational data is taken from Marconi et al. (1996). 
The set of lines is not as complete but is probably sufficient 
for estimating Uo- Marconi et al. reported that the excita- 
tion was lower than in Circinus. A similar result is found 
here. A value four times lower than in Circinus leads to a 
reasonable fit to the gas excitation as revealed by the S and 
Si line ratios. The fit is shown in Fig. lb (Model B: U =0.125 
and /^d=0.1). Systematic but equal departures from unity 
(e.g. the S lines) might partly be due to peculiar abun- 
dances relative to solar. In this model, we arbitrarily adopted 
the same density as for Circinus since no published mea- 
surements of the density from the [Nev] 14.3 /^m/24.3 /im 
ratio was available at the time of writing. The excitation 
indicated by the ratio [Fevn]/[Fex] is not so badly fitted 
but [Fexi] A7892 is certainly discrepant. As was the case for 
Circinus, assigning more than 50% of the radiation pressure 
to cloud acceleration (i.e. >J/o>0.5) was found to worsen 
the fit. On the other hand, the amount of dust allowed did 
not appear to strongly affect the fit as long as /id< 0.5. A 
good fit to the value of a' as a function of the dust content 
is a' = (jio/1200)(Uo/0.125) [3.1 x 10~ 7 + 5.8 x lO^Vu] ■ 

If we vary no and keep the ionization parameter con- 
stant, by equation (1) qo must scale with no. Therefore the 
modified accelerations a' will also scale linearly with no since 
radiation pressure scales with ionizing flux. Note that a sim- 
ilar model to that of Fig. lb but for which the dust content 
has been increased to ud=0.3 (i.e. Model C in Table 1), the 
acceleration a' becomes comparable to that of Model A with 
Uo four times larger (but three times less dust). 

For Models A and B, the contrast between the density 
at the irradiated face and the 'emissivity averaged' value, 
n/no, is 3.75 and 1.75 respectively. Table 1 list other relevant 
quantities for the two models shown in Fig. 1 as well as for 
a dustier Model C characterized by ud=0.3. 



4 RESULTS 

For an empirically determined value of Uo, a 1 will scale lin- 
early with the ionizing flux qo. Therefore, for a given value of 
Uo , the predicted value of a 1 scales linearly with the assumed 
value of the density no. Using the [Nev] 14.3 /im/24.3 (im 
line ratio, one can determine the density of the coronal line 
region and therefore a'. It is found that the mean density 
no as defined in Paper I for a pressure stratified slab is the 
same as that inferred from the calculated far-infrared [Ne v] 
line ratio. Therefore this line ratio reveals us the 'charac- 
teristic' density of the highly excited component. Having 
determined empirically the ionization parameter character- 
izing the coronal gas in two AGN, we can now proceed to 
compare the calculated a' of models with, g pu ii, the strength 
of the gravitational field expected in a galactic nucleus. 
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Figure 2. A plot of the gravitational pull, g pu ii, exerted by the 
nuclear mass distribution of the Milky Way (dash line) as a func- 
tion of galactic radius. The continuous and the dash-dotted lines 
consider the effect of a m ore massive BH of 10 7 or 3 x 10 7 M Q , 
respectively (see Section fl l[ ). The modified acceleration a' due 
to radiation pressure for clouds of different densities n is indicated 
using horizontal lines (whatever n, it is assumed that Uo=0.5). 
The minimum radius for which coronal gas emission can be ac- 
celerated (i.e. a' > g pu ii) is indicated in the right margin for the 
case of a 10 7 Mq BH. 



4.1 The nuclear mass model: M t 



+ bh 



(r) 



For radiative acceleration to be responsible of the blueshift 
of the coronal lines, it must overcome the gravitational pull 
of the combined mass of the central stars and the nuclear 
BH. To illustrate the force entailed, let us adopt the bet- 
ter studied Milky Way mass distribution, M t+ i,h(r), inside 
a radius r < 200 pc. From 0.1 to 2pc we borrow the mass 
model of Genzel et al. (1997) which we join at 2pc to the 
isotropic orbit stellar mass model of Lindqvist et al. (1992; 
cf. their Table 3) which is based on the projected mass 
method for a dominant central mass. The resultant grav- 
itational pull is simply g pu u = GM ir+ th/i r2 and is plotted 
in Fig. ^ (dash line). One interpretation proposed by Gen- 
zel et al. for the innermost mass is that it might consist of 
a BH of w 2.8 x 10 6 A'/ Q . Note that beyond r £ 30 pc the 
gravitational pull from the stars dominate that of the BH. 
On account of the much stronger nuclear activity character- 
izing Seyferts, it is plausible that a typical AGN could have 
a more massive BH than the Milky Way. For instance, esti- 
mates of the central mass in NGC 1068 ranges from 10 7 Af n 
(iGrccnhill fe Gwinn 1997| ) to 3 x lO 7 Af ( |Gallimore et al 



1996). In an attempt to allow for this possibility, the effect 



of simply substituting a more massive point source is shown 
in Fig. H in which the continuous and the dotted-dash curves 
correspond to BH masses of 10 7 Mq and 4 x 10 7 Mq, respec- 
tively. 
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Table 1. Parameters characterizing the models. 



Model 


U 


no 


Md 


JV t /10 21 


n 


Ay 


511 


A 


0.50 


1200 


0.10 


6.2 


4500 


0.15 


1.22 


B 


0.125 


1200 


0.10 


3.0 


2100 


0.07 


0.51 


C 


0.125 


1200 


0.30 


2.6 


2000 


0.19 


0.25 



//lO" 6 


n/no 






P/fc/lO 7 ^ 


L (pc) 


2.63 


3.75 


0.06 


-0.36 


13.2 


0.68 


0.89 


1.75 


0.22 


-0.46 


7.8 


0.55 


2.05 


1.65 


0.04 


-0.18 


8.1 


0.35 



a Fraction of ionizing photons escaping the slab unabsorbed. 

b Ratio of dynamical acceleration to total radiative acceleration (Mathews 1986). 

c Isotropic pressure of the ambient cloud-confining medium. It would scale with no if Uo is kept constant. 
d Geometrical depth of the slab in parsecs. L would scale inversely with no if Uo is kept constant. 



4.2 Minimum NLR radius 

It is not possible at this point to realistically assess how 
different the nuclear stellar field would be in either more 
massive galaxies than ours or in nuclei with more massive 
BHs. With this caveat in mind, let us adopt the continuous 
curve in Fig. ^ of a 10 7 Mq BH as being representative of low 
luminosity AGN. After determining the different values of a' 
expected if we simply vary no of Model A (with Uo=0.5), we 
can proceed to derive the minimal radius beyond which the 
clouds have their acceleration a' exceed g pu ii. The allowed 
domain in r for uniform acceleration of the coronal gas lies 
to the left of the g pu ii curve in Fig. |^ as shown by the hor- 
izontal solid lines, each corresponding to different values of 
the mean densities n. For the particular case of the Circinus 
galaxy for which we know the density, a minimum radius of 
12 pc is inferred inside of which the coronal gas would fall 
or stall instead of being radially accelerated. The true posi- 
tion of the nucleus in Circinus is not known with precision 
since its nucleus is heavily reddened and the true kinematic 
centre may not correspond to the continuum peak. We can- 
not yet therefore verify the validity of r m j n . However, in the 
case of NGC 1068, Marconi et al. (1996) found that the in- 
frared lines were displaced by 0.5" from the optical contin- 
uum peak, which corresponds to a nuclear offset of ~ 40 pc 
(D=18Mpc). The gas density is not yet determined but if 
comparable to Circinus, we would expect a much larger r m i n 
in NGC 1068 based on its lower excitation by a factor four. 

The main conclusions are the following. Provided there 
is dust left out in the gas, radiation pressure not only can 
induce a density gradient in the photoionized clouds but 
also be responsible for their outward acceleration despite the 
presence of strong gravitational field. It is found that higher 
values than ud=0.1 are detrimental to the fit in the case of 
Circinus. Based on empirically measured densities and ex- 
citation, a minimum nuclear distance r m i n is predicted to 
characterize the observable coronal gas. The higher the den- 
sity, the smaller r m i n . The advent of high spatial resolution 
line mapping will allow testing whether such a correspon- 
dence between r and n exists. Finally, because of the large 
value of Uo, the fraction of ionizing photons escaping the 
clouds can be quite small even with small amounts of dust. 
By empirically setting the thickness of the clouds to be such 
that the observed [Si vi]/[Siix] ratio is reproduced by the 
model, one finds that only 6% of the impinging photons in 
Model A are not absorbed (instead of 35-40% in the dust- 
free case). One may therefore conceive the possibility that 
the coronal clouds are not matter bounded but marginally 
thick and bounded by dust absorption. The most inner nu- 
clear regions (presumably of higher visibility in Seyferts 1) 



which still emit coronal lines may correspond to the missing 
domain between the BLR and the NLR (this gas would have 
to be proportionally denser to have a 1 > g pu a)- In effect, as 
shown by Netzer & Laor (1993), dust embedded photoion- 
ized gas of high excitation presents a greatly reduced line 
efficiency explaining the relative line weakness of this puta- 
tive intermediate region. 
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